Patent Application No. 10/697,362 

Amdt. Dated August 7, 2007 

Response to Advisory Action of July 16, 2007 



Docket No. SAM2.PAU.33 



Amendments to the Claims: 

1.-4. (Canceled). 

5. (Previously presented) A method for reducing noise in a digital image formed 
from a plurality of pixels including a given pixel, the method comprising the steps of: 

computing global statistics from the image; 

computing local statistics for the given pixel; 

configuring a local filter using the local and global statistics; 

filtering the given pixel using the local filter to reduce image noise; 

wherein the step of computing the local statistics for the given pixel 
further includes the steps of: 

selecting a window containing the given pixel and a plurality of 
neighboring pixels; 

computing a 2-D local variance of the given pixel based on information 
related to the pixels in the window; 

computing a plurality of 1-D local variances along multiple directions 
through the given pixel in the window; and 

detecting a local edge direction by selecting one of the directions with the 
smallest 1-D local variance, 

wherein the step of configuring the local filter using the local and global 
statistics includes the steps of: 
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selecting the detected local edge direction L as the direction of the local 

filter; 

for the detected local edge direction L, computing a 1-D filter strength as a 
function of a square root of a local variance and a global noise standard deviation; 

computing a 2-D filter strength as a function of the local variance and the 
global noise standard deviation; and 

configuring the local filter for the detected local edge direction L based on 
the 1-D and 2-D filter strengths. 



6. (Previously presented) The method of claim 5, wherein the step of configuring 
the local filter using the local and global statistics further includes the steps of: 

selecting the detected local edge direction Lk (k = 1, 2, 3, or 4) as the 
direction of the local filter; 

for the detected local edge direction L k , computing the 1-D filter strength 
a k = min(2<T,max(3cr - a k fi))/(2a) wherein a is the global noise standard deviation, 

and a k is the square root of the local variance (k = 1, 2, 3, or 4); 

computing the 2-D filter strength a 0 = min(2cr, max(3cr - cr 0 ,0))/(2cr) 
wherein a 0 is a square root of a 2-D local variance; and 

configuring the local filter f k for the detected local edge direction Lk 

according to the following conditions, wherein a k is the filter strength along the local 
edge direction L k [k = 0 (non-edge), 1 (horizontal), 2 (vertical), 3 (upper left to lower 
right, 4 (upper right to lower left)]: 
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(i) if the detected direction is Lj, then /, is configured as a 2-D local 



filter for horizontal direction, wherein: 



"5 



a 0 +3a,(l-a 0 ) a 0 +3(3-2a,)(l-« 0 ) a 0 + 3a,(l-a 0 ) 



(ii) if the detected direction is L2, then f 2 is configured as a 2-D local 



filter for vertical direction, wherein: 



a 0 a 0 + 3a 2 (l-a 0 ) a 0 
a 0 a 0 + 3(3-2a 2 )(l-a 0 ) a 0 
a 0 a 0 +3a 2 (l-a 0 ) a 0 



(iii) if the detected direction is L3, then / 3 is configured as a 2-D local 
filter for the diagonal direction from upper left to lower right, wherein: 



h 9 



a 0 + 3a 3 (l-a 0 ) 



a 0 + 3(3-2a 3 )(l-a 0 ) 



a 0 +3a 3 (l-a 0 ) 



; and 



(iv) if the detected direction is L4, then f 4 is configured as a 2-D local 
filter for the diagonal direction from upper right to lower left, wherein: 



h 9 



a 0 +3a 4 (l-a 0 )' 



a 0 +3(3-2a 4 )(l-a 0 ) 



a 0 +3a 4 (l-a 0 ) 
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7. (Previously presented) A method for reducing noise in a digital image formed 

from a plurality of pixels including a given pixel, the method comprising the steps of: 

computing global statistics from the image; 

computing local statistics for the given pixel; 

configuring a local filter using the local and 
global statistics; 

filtering the given pixel using the local filter to reduce image noise; 

wherein the step of computing the local statistics for the given pixel 
further includes the steps of: 

selecting a window containing the given pixel and a plurality of 
neighboring pixels; 

computing a 2-D local variance of the given pixel based on information 
related to the pixels in the window; 

computing a plurality of 1-D local variances along multiple directions 
through the given pixel in the window; and 

detecting a local edge direction by selecting one of the directions with the 
smallest 1-D local variance, 

wherein the step of computing the global statistics comprises the step of 
estimating a global noise standard deviation 6 to generate the global statistics, 

wherein the step of estimating the global noise standard deviation a 
further includes the steps of: 

dividing the input image into overlapping or non-overlapping blocks; 

computing a mean and a standard deviation for each block; 
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finding the smallest standard deviation d 0 and its corresponding mean 
detecting block saturation due to noise; 

compensating for the smallest standard deviation d 0 to generate a 
compensated smallest standard deviation d 0 ; 

selecting the block standard deviations d n that are within a range of the 
compensated smallest standard deviation d 0 ; and 

averaging the selected block standard deviations d n to generate an 
estimate of the global noise standard deviation o . 

8. (Previously presented) The method of claim 7, wherein the block size is 7 x 7 
or 5x9 pixels. 

9. (Previously presented) The method of claim 7, wherein the steps of detecting 
block saturation due to noise and compensating for the smallest standard deviation 
include the steps of: 

defining an upper pixel value limit UL, a lower pixel value limit LL, and a 
mid value M between UL and LL, 

wherein if the mean m 0 is less than the mid value M, and the smallest 

standard deviation is greater than a difference between the mean m 0 and the lower limit 
LL 9 then saturation has occurred at the lower limit LL, and the smallest standard deviation 
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d Q is compensated by adding thereto a compensation term that is a function of the 
smallest standard deviation d 0 and said difference between the mean m 0 and the lower 

limit LL, to generate the compensated smallest standard deviation d 0 ; 

else if the mid value Mis less than the mean m 0 , and the smallest standard 
deviation d Q is greater than the difference between the upper limit UL and the mean m 0 , 
then saturation has occurred at the upper limit UL, and the smallest standard deviation d Q 
is compensated by adding thereto a compensation term that is a function of the smallest 
standard deviation d 0 and the difference between the upper limit UL and the mean m 0 , to 

generate the compensated smallest standard deviation d 0 ; 

otherwise, no saturation has occurred, wherein d 0 =d Q . 

10. (Previously presented) The method of claim 7, wherein the steps of detecting 
block saturation due to noise and compensating for the smallest standard deviation 
include the steps of determining the following: 

where UL is an upper pixel value limit, LL is a lower pixel value limit, M 
is a mid pixel value, and UL <M<LL, if the mean m 0 <M and the smallest standard 
deviation d 0 >m 0 -LL , then saturation has occurred at the lower limit LL, wherein d Q is 
compensated as d 0 = d 0 + K • (d 0 - (m 0 - LL)) , such that K is a compensation factor; 



7 of 16 



Patent Application No. 10/697,362 Docket No. SAM2.PAU.33 

Amdt. Dated August 7, 2007 

Response to Advisory Action of July 16, 2007 

else if the mean m 0 >M and the smallest standard deviation 
d Q >UL-m 0 , then saturation has occurred at the upper limit UL, wherein d 0 is 
compensated as d 0 =d 0 +K- (d Q - (UL - m 0 )) ; 

otherwise, no saturation has occurred, wherein d 0 =d Q . 

11. (Previously presented) The method of claim 10, wherein LL = 0,UL = 255, 
and M= 128. 

12. (Previously presented) The method of claim 7, wherein the step of selecting 
the block standard deviations includes the steps of selecting the block standard deviation 
d n for averaging if | d n - d 0 |< max((i 0 ,l) . 



13.-16. (Canceled). 



17. (Previously presented) A noise reduction system for reducing noise in a digital 
image comprising pixels, the system comprising: 

a global statistics module that computes global statistics from the image; 
a local statistics module that computes local statistics for each of a 
plurality of image pixels including a given pixel; 

a filter configuration module that uses the local and global statistics for 
the given pixel to configure a local filter for filtering the given pixel; 
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the local filter as configured by the filter configuration module, being 
adapted for filtering the given pixel to reduce image noise, wherein the local statistics 
module computes the local statistics for the given pixel by: 

selecting a window containing the given pixel and a plurality of 
neighboring pixels; 

computing a 2-D local variance of said pixel based on information related 
to the pixels in the window; 

computing a plurality of 1-D local variances along multiple directions 
each defined by a pair of the pixels in the window; and 

detecting a local edge direction for the given pixel by selecting one of the 
directions with the smallest 1-D local variance, 

wherein the global statistics module estimates a global noise standard 
deviation a to generate the global statistics, 

wherein the filter configuration module configures the local filter for each 
pixel using the local and global statistics by: 

selecting a detected local edge direction L as the direction of the local 

filter; 

for the detected local edge direction L computing a 1-D filter strength as a 
function of the square root of the local variance and the global noise standard deviation; 

computing a 2-D filter strength as a function of a local variance and a 
global noise standard deviation; and 

configuring the local filter for the detected local edge direction L based on 

the 1-D and 2-D filter strengths. 
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18. (Previously presented) The system of claim 17, wherein the filter 
configuration module configures the local filter for each pixel using the local and global 
statistics by: 

selecting the detected local edge direction Lj< (k = 1, 2, 3, or 4) as the 
direction of the local filter; 

for the detected local edge direction Lk , 

computing the 1-D filter strength a k - min(2o*, max(3cr - a k ,0)) /(2a) 
wherein a is the global noise standard deviation, and a k is the square root of the local 

variance (k = 1, 2, 3, or 4), and the 2-D filter strength 

a 0 = min(2cr, max(3<j - ct 0 ,0))/(2<j) wherein a 0 is a square root of a 2-D 

local variance; and 

configuring the local filter f k for the detected local edge direction L k 

according to the following conditions: 

(i) ifthe detected direction is Li, then /, is configured as a 2-D local 
filter for horizontal direction, wherein: 

a 0 a 0 
ff 0 +3(3-2a,)(l-a 0 ) a 0 +3a,(l-a 0 ) ; 
a 0 a 0 

(ii) if the detected direction is L 2 , then f 2 is configured as a 2-D local 
filter for vertical direction, wherein: 



Jl 9 



a 0 +3a 1 (l-or 0 ) 
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f= l 

9 



a 0 a 0 +3a 2 (l-a 0 ) 



« 0 +3(3-2a 2 )(l-a 0 ) a 0 



_a 0 a 0 +3a 2 (l-a 0 ) 



a, 



o _ 



(iii) if the detected direction is L3, then f 3 is configured as a 2-D local 
filter for the diagonal direction from upper left to lower right, wherein: 



a 0 +3a 3 (l-« 0 ) 



a 0 + 3(3-2a 3 )(l-a 0 ) 



a 0 +3a 3 (l-« 0 )_ 



; and 



(iv) if the detected direction is L4, then f 4 is configured as a 2-D local 
filter for the diagonal direction from upper right to lower left, wherein: 



9 



a 0 +3a 4 (l-a 0 )" 



a 0 +3(3-2a 4 )(l-a 0 ) 



a 0 +3a 4 (l-« 0 ) 



19. (Previously presented) A noise reduction system for reducing noise in a digital 
image comprising pixels, the system comprising: 

a global statistics module that computes global statistics from the image; 

a local statistics module that computes local statistics for each of a 
plurality of image pixels including a given pixel; 

a filter configuration module that uses the local and global statistics for 
the given pixel to configure a local filter for filtering the given pixel; 

the local filter as configured by the filter configuration module, being 
adapted for filtering the given pixel to reduce image noise, wherein the local statistics 
module computes the local statistics for the given pixel by: 
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selecting a window containing the given pixel and a plurality of 
neighboring pixels; 

computing a 2-D local variance of said pixel based on information related 
to the pixels in the window; 

computing a plurality of 1-D local variances along multiple directions 
each defined by a pair of the pixels in the window; and 

detecting a local edge direction for the given pixel by selecting one of the 
directions with the smallest 1-D local variance, 

wherein the global statistics module estimates a global noise standard 
deviation a to generate the global statistics, 

wherein the input image comprises a plurality of blocks, and wherein the 
global statistics module further comprises: 

a mean and standard deviation module that computes the mean and the 
standard deviation for each block ; 

a minimum finder module that finds the smallest standard deviation d 0 and 

its corresponding meanm 0 ; 

a saturation detector that detects block saturation due to noise; 

a saturation compensator that compensates for the smallest standard 

deviation^ to generate a compensated smallest standard deviation d 0 ; and 

a selective averaging module that selects the block standard deviations 

d n that are within a range of the compensated smallest standard deviation d 0 , and 
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averages the selected block standard deviations d n to generate an estimate of the global 
noise standard deviation a . 

20. (Previously presented) The system of claim 19, wherein the block size is 7 x 7 
or 5 x 9 pixels. 

21. (Previously presented) The system of claim 19, wherein: 

an upper pixel value limit is denoted UL, a lower pixel value limit is 
denoted LL, and a mid value M is between UL and LL, 

wherein the saturation detector determines if the mean m 0 is less than the 

mid value M 9 and the smallest standard deviation is greater than a difference between the 
mean m 0 and the lower limit LL, indicating that saturation has occurred at the lower limit 
LL 9 and if so, the saturation compensator compensates for the smallest standard deviation 
d 0 by adding thereto a compensation term that is a function of the smallest standard 

deviation d Q and said difference between the mean m 0 and the lower limit LL, to 

generate the compensated smallest standard deviation d 0 ; 

else if the saturation detector determines that the mid value M is less than 
the mean m 0 , and the smallest standard deviation d Q is greater than a difference between 

the upper limit UL and the mean m 0 , indicating saturation has occurred at the upper limit 
UL, the saturation compensator compensates for the smallest standard deviation d 0 by 
adding thereto a compensation term that is a function of the smallest standard deviation 



13 of 16 



Patent Application No. 10/697,362 

Amdt. Dated August 7, 2007 

Response to Advisory Action of July 16, 2007 



Docket No. SAM2.PAU.33 



d Q and said difference between the upper limit UL and the mean m Q , to generate the 
compensated smallest standard deviation d Q ; 

otherwise, no saturation has occurred, wherein d 0 =d 0 . 

22. (Previously presented) The system of claim 21, wherein: 

if the mean m 0 <M and the smallest standard deviation d 0 >m 0 -LL , 

indicating saturation has occurred at the lower limit LL, then d 0 is compensated for as 

d 0 =d 0 +K -(d 0 - (m 0 - LL)) , such that K is a compensation factor; 

else if the mean m 0 >M and the smallest standard deviation 
d 0 >UL-m Q9 indicating saturation has occurred at the upper limit UL, then d 0 is 
compensated for as d 0 =d 0 +K- (d 0 - (UL - m 0 )) ; 

otherwise, no saturation has occurred, wherein d 0 =d Q . 

23. (Previously presented) The system of claim 22, wherein LL = 0, UL = 255, 
and M= 128. 

24. (Previously presented) The system of claim 19, wherein the block standard 
deviations d n are selected for averaging if | d n - d 0 \< max(d 0 ,l) . 
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